Pair Creation and an X-ray Free Electron Laser 
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Using a quantum kinetic equation coupled to Maxwell's equation we study the possibility that 
focused beams at proposed X-ray free electron laser facilities can generate electric field strengths 
large enough to cause spontaneous electron-positron pair production from the QED vacuum. 
Our approach yields the time and momentum dependence of the single particle distribution 
function. Under conditions reckoned achievable at planned facilities, repeated cycles of particle 
creation and annihilation take place in tune with the laser frequency. However, the peak particle 
number density is insensitive to this frequency and one can anticipate the production of a few 
hundred particle pairs per laser period. Field-current feedback and quantum statistical effects 
are small and can be neglected in this application of non-equilibrium quantum mean field theory. 
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The QED vacuum is unstable in the presence of a 
strong external field and decays by emitting electron- 
positron pairs. The pair production rate was first calcu- 
lated for a static, homogeneous electric field in the early 
part of the last century |^ and since then many aspects 
have been studied in detail. Using the Schwinger for- 
mula Ql one finds that a sizeable rate requires a field: 
Ecr ■— ml/e = 1.3 x 10^^ V/m, which is very difficult to 
produce in the laboratory. (We use ?i = c = 1.) 

The possibility of spontaneous pair creation from the 
vacuum is of particular interest in ultra-relativistic heavy 
ion collisions ||, ^ . Since the QCD string tension is large 
{y/a ~ 400 MeV), fiux tube models of the collision gen- 
erate a background field that is easily strong enough to 
initiate the production process via a Schwinger-like mech- 
anism. Feedback between the external field and the field 
created by the produced particles' motion drives plasma 
oscillations. This is the back-reaction phenomenon, 
which has been much discussed ||, |[ ^. While this 
Vlasov-equation-based approach has met with some phe- 
nomenological success, a rigorous justification in QCD is 
wanting. 

Pair creation by laser beams in QED has also been 
discussed [Q, ^, and proposed X-ray free electron laser 
(XFEL) facilities at SLAC § and DESY |o), which 
could generate field strengths |l^ « O.lEcr, 
promise to provide a means to explore this phenomenon. 

Vacuum decay is a far-from-equilibrium, time-depen- 
dent process and hence kinetic theory provides an ap- 
propriate descriptive framework. For spatially homoge- 
neous fields, a rigorous connection between kinetic the- 
ory and a mean-field treatment of QED has been es- 
tablished ijl^, The derivation makes plain that the 
true kinetic equation's source term is intrinsically non- 
Markovian, and this is expressed in properties of the so- 
lution im Herein we use this quantum Vlasov 



equation to obtain a description of the time evolution of 
the momentum distribution function for particles pro- 
duced via vacuum decay at the planned XFEL facilities. 

A gauge and Lorentz invariant description of an elec- 
tromagnetic (e.m.) field is obtained using 



J" = -F^^Ff"" = -(E^- 
4 ^ 2^ 



g = ^F^.F^-' = E 



B. 

(1) 



An e.m. plane wave always fulfills T — Q = Q and such 
a light-like field cannot produce pairs [|7|. Therefore, to 
produce pairs it is necessary to focus at least two coher- 
ent laser beams and form a standing wave. Subsequently 
we assume an "ideal experiment:" owing to the diffrac- 
tion limit the spot radius of the crossing beams cannot be 
smaller than the wavelength, so we choose r„ « A; and 
we assume a space volume in which the electric field is 
nonzero but the magnetic field vanishes. Even for care- 
fully chosen X-ray optical elements, such a situation is 
impossible to achieve in practice, which means that the 
field strength actually available to produce particles is 
weaker than the peak field value and hence our estimate 
of the production rate will be an upper bound. 

Our idealised experiment is realised through a vec- 
tor potential A^{t) that generates an antiparallel elec- 
tric field Eext(t) = —dA{t)/dt: this is the background 
field, which in temporal gauge, Aq — 0, reads A{t) = 
{0,0, A{t)), and the alternating laser field is then 



Ee,t{t)^ (0,0, Eosinint)); B = 0, 



(2) 



where fl = 27r/A is the frequency. (Table | provides our 
laser field parameters.) The total electric field is E{t) = 
Eext{t) + Eint{t) where the internal field is induced by 
the back-reaction mechanism. As we will show, Eint can 
be neglected under anticipated XFEL conditions. 
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TABLE I: Laser field parameters are specified in columns one 
and two: Set I is XFEL-like; Set II is strong. Columns three 
and four describe the density, nmax{t>), and total number 
of produced particles, N{t^) — X^Umaxity), Eq. (^), where 
is the time at which the number density reaches its (for 
weak fields, local) maximum. A typical laser pulse length is 
~ 80 fs 0. 





A (nm) 


Eo (V/m) 




(fnr^) 


N(t>) 


Set la 


0.15 


1.3 X 10'' 


4.6 X 10" 


13 


~ lO'' 


Set lb 


0.075 


1.3 X 10'^ 


4.6 X 10" 


13 


~ 10^ 


Set Ila 


0.15 


1.3 X 10'® 


7.2 X 10" 


10 


~ 10'^ 


Set lib 


0.075 


1.3 X 10'® 


6.4 X 10" 


10 


~ 10^ 



A key quantity in the description of nonequilibrium 
particle production processes is the single-particle mo- 
mentum distribution function, f{p, t). For Dirac particles 
coupled to an Abelian gauge field it can be obtained by 
solving a quantum Vlasov equation ||l3|, Once f{p, t) 
is known the calculation of the produced particle number 
density and total particle number is straightforward. In 
general the Vlasov equation involves source and collision 
terms, and its coupling to Maxwell's equation provides 
for the field-current feedback typical of plasmas. The 
absolute and relative importance of these terms depends 
on the magnitude of the background field and the mass 
of the produced particles §. For the relatively weak 
XFEL-like fields we expect the produced-particle num- 
ber density to be small and hence collisions to be rare. 
Therefore we neglect the collision term and arrive at the 
following pair of coupled equations: 



df{p,t) 
dt 



,E{t)el /•*^^,ei?(<')[l-2/(p,t')] 



X cos 



dr uj{p, t) 



(3) 



u;^p) df{p,t) eE{t)el 



zE{t)p\\ dt 8cj4(p)p|| 



(4) 



with the three- vector momentum p= the trans- 

verse mass-squared e\ = rnl+p\; and the total energy- 



squared Uj'^{p,t) = £2 +p2. 



P3 



'.A{t). 



The effect of quantum statistics on the particle pro- 
duction rate is evident in the factor "[1 — 2 /]" in Eq. (||), 
which ensures that no momentum state has more than 
one spin-up and one spin-down fermion. In addition, 
both this statistical factor and the "cos" term introduce 
non-Markovian character to the system: the first couples 
in the time history of the distribution function's evolu- 
tion; the second, that of the field. One anticipates that 
such memory effects are only important for very strong 
fields: E ~ Ecr] i.e., when the time-scale characteristic 
of the produced particles' Compton wavelength is simi- 
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FIG. 1: Time and momentum dependence of the distribu- 
tion function /(p, t) for p_i_ — obtained with the conditions 
described in Set la of Table ^ This alternating laser field pro- 
duces repeated cycles of particle creation and annihilation. 



lar in magnitude to the time taken to tunnel through the 
barrier [l^. We will verify that such effects are 

unimportant for XFEL-like fields and that, in this ap- 
plication, one may employ the low density (l.d., / <C 1) 
approximation to distribution function: 



to 



'dr^-^ rdt'4^ 

2cj^p,T)Jt^ LO^{p,t') 
t 

2 / dr' uj{p,t') 



X cos 



(5) 



Equations (g) and are coupled: Eq. (3) depends 
on the total field, E{t), calculated from Eq. (4), and the 
electric field itself depends on the polarisation and con- 
duction currents, which are proportional to df /dt and 
/, respectively. For strong fields the internal currents 
are non-negligible and field-current feedback will gener- 
ate plasma oscillations. However, as we will demonstrate, 
that is not relevant for realisable XFEL fields. 

We turn now to a quantitative analysis of two exem- 
plary electric field strengths in Eq. (||), solving Eqs. (j^) 
and (|) using a 4 order Runge-Kutta method. The 
weaker field: E = 0.1 Ecr', i.e.. Sets I in Table |, should 
be obtainable at the proposed XFEL facilities lOj. 
Sets II in the table, obtained with E = E^r, provide a 
strong field comparison. 

Figure ^ depicts the time evolution of the momentum 
distribution function for electrons produced under cur- 
rently anticipated XFEL conditions. This relatively weak 
alternating field generates repeated cycles of pair cre- 
ation, with about 1.6 x 10^ local maxima in one laser 
pulse. The production rate: S(j),t) = df(p,t)/dt, is 
clearly time-dependent and vanishes at these local max- 
ima. Therefore an estimation of S' at a fixed time, as 
in Refs. |l^, |l^, is potentially misleading. Under these 
conditions there is no accumulation of produced particles 
because creation and annihilation processes are balanced; 
i.e., all the particles produced have sufficient time to an- 
nihilate again before the next cycle of particle creation. 
Hence the number of particles cannot exceed the value 
at the first of the repeated maxima. 
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E„=Ec, X,=0.15nm 
E„=E;. , X=0.075nm 
E„=0.1Ec,A.=0.15nm 
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FIG. 2; Time evolution of the number density for Sets I 
and II, Table |. In strong fields, particles accumulate, leading 
to the almost complete occupation of available momentum 
states. In weak fields, repeated cycles of particle creation and 
annihilation occur in tune with the laser frequency: flia = 
2.0 X 10^** 2txs-^; fiib = 4.0 x 10^** 2txs-^. 



The situation is markedly different for strong fields, as 
made plain by Fig. |^, which depicts the time evolution 
of the number density: 



n{t) 



(2^ 



(6) 



When E = Ecr the interval between bursts of particle 
creation is so small that not all the particles can anni- 
hilate before the next cycle of particle production be- 
gins. Therefore the particles accumulate until almost all 
momentum states are completely occupied. In this case 
field-current feedback effects may become important. 

Having calculated the time-dependent number density 
we can estimate the number of particles produced in the 
spot volume via vacuum decay. Assuming the diffraction 
limit on focusing can be reached then = A, Va ~ X^, 
and the peak particle number is 



7V(i>)«A3n(t>), 



(7) 



where t> is the time-location of the peak number den- 
sity (which is any of the local maxima for XFEL fields). 
Our results are presented in Table |. We remark that, 
following a very different route, Ref. also arrives at 
a value similar to our Set la estimate of ~ 10'^ particles. 
Increasing the peak field strength by one order of mag- 
nitude increases the produced particle yield by 10'^, as is 
obvious from Fig. ^. 

Our complete numerical solution also provides the full 
momentum distribution, which is depicted in Figs. ^. For 
both weak and strong fields the maximum transverse mo- 
mentum is approximately 1 MeV but the maximum longi- 
tudinal momentum is much larger for strong fields. This 
result owes itself to the minimal coupling of the vector 
potential to the parallel momentum component but is 
invisible for weak fields. Furthermore the presence of 
non-Markovian effects for strong fields (see discussion of 
Eq. (H)) is evident in the irregular structure in the lower 
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FIG. 3: Momentum-dependence of the distribution function: 
Upper panel. Set la; Lower panel Set Ila, Table |. Note the 
vastly different magnitudes. The accessible range of longitu- 
dinal momenta is magnified for strong fields and the irregular 
structure expresses the influence of non-Markovian effects in 
the source term for this case. 



panel cf. the smooth distribution in the XFEL-like field 
distribution depicted in the upper panel. 

In Fig. U we plot the time evolution of the current. 
For an XFEL-likc field the induced current is zero on the 
scale of this figure. (NB. It is a factor of 10^ smaller than 
the external current.) The plasma oscillation period in 
this system can be estimated using the ultrarelativistic 
formula 



IJR 



^2 



.e/2 



6.5 X lO^i^s 



(8) 



cf. the laser periods: 2.5 x 10~^^ s and 5.0x10"^^ s, which 
shows that the oscillatory behaviour evident in the fig- 
ure is prescribed by the laser frequency. (NB. It em- 
phasises, however, that were the amplitude of the in- 
duced current larger, field-current feedback may become 
important.) The behaviour of the strong field current 
is instructive. The focused laser beams generate a vol- 
ume of high field strength in which particles are created. 
The strong field persists long enough for all the parti- 
cles to be accelerated to their maximum velocity, almost 
the speed of light in this case, which explains the ap- 
pearance of the plateaux: the current saturates when all 
the produced particles reach their terminal velocity. The 
plateaux end as the alternating laser field changes sign 
and accelerates the electrons in the opposite direction. 
Each successive plateau is higher because, as explained 
in connection with Fig. ^, for strong fields the period be- 
tween particle production cycles is small and hence there 
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FIG. 4: Current for Sets II of Table ^ The plateaux occur 
when the produced particles reach their terminal velocity, al- 
most the speed of light. The oscillatory behaviour is driven by 
the laser frequency. (NB. For an XFEL-like field the current 
is zero on the scale of this figure.) 



is a net accumulation of particles. 

We have explored the possibility of electron-positron 
pair production through vacuum decay at the focus of 
achievable X-ray free electron laser beams via the self- 
consistent solution of a quantum Vlasov equation cou- 
pled to Maxwell's equation. This is a rigorous appli- 
cation of non-equilibrium quantum mean field theory. 
Our complete calculation confirms that the effect will 
be observable if proposed facilities achieve their design 
goals. Under such conditions the production rate is time- 
dependent, with repeated cycles of particle production 
and annihilation in tune with the laser frequency, but 
the peak particle number is independent of the laser fre- 
quency: up to 10^ pairs may be produced in the spot 
volume. Our procedure yields the full single particle mo- 
mentum distribution function (Fig. ^ and hence provides 
comprehensive phase space information about the e+e^ 
pairs. It also verifies that, at present design goals: col- 
lisions among the produced particles can be neglected 
because the maximum particle number density is low 
(Fig. 1^); a Schwinger-like source term is adequate because 
the non-Markovian features of the true source term have 
no impact (Fig. |^) and hence the low density approxi- 
mation to the distribution function, Eq. (||), is accurate; 
and field-current feedback will not influence the produc- 
tion process (Fig. 0). 
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